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Abstract. We discuss the angular correlation function of Sunyaev-Zel'dovich (SZ)-detected galaxy clusters as a 
cosmological probe. As a projection of the real-space cluster correlation function, the angular function samples the 
underlying SZ catalog redshift distribution. It offers a way to study cosmology and cluster evolution directly with 
the two-dimensional catalog, even before extensive follow-up observations, thereby facilitating the immediate 
scientific return from SZ surveys. As a simple illustration of the information content of the angular function, 
we examine its dependence on the parameter pair (QmjOs) in fiat cosmologies. We discuss sources of modeling 
uncertainty and consider application to the future Planck SZ catalog, showing how these two parameters and the 
normalization of the SZ flux-mass relation can be simultaneously found when the local X-ray cluster abundance 
constraint is included. 
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1. Introduction 

The Sunyaev-Zel'dovich (SZ) effect (Sunyaev & Zel'dovich 
1970, 1972) has become a practical observational tool 
for studying galaxy clusters and cosmology (for recent 
reviews see Birkinshaw 1999; Carlstrom et al. 2002). 
Current observations of individual clusters, when com- 
bined with X-ray observations, constrain cosmological pa- 
rameters via gas mass fractions (Grego et al. 2002) and 
angular-diameter distance determinations (Mason et al. 
2001; Jones et al. 2001; Reese et al. 2002). Multi-band mil- 
limeter observations of a handful of clusters have already 
been used to set limits on peculiar velocities (Holzapfel 
et al. 1997; Benson et al. 2003), and theoretical studies of 
this technique show its promise for the future (Aghanim et 
al. 2002a, 2002b; Holder 2002). A new generation of opti- 
mized, dedicated instruments, both large-format bolome- 
ter arrays and interferometers with high sensitivity re- 
ceivers, will qualitatively improve these studies. And the 
arrival of these instruments within the next few years, in 
addition to the Planck mission^, will move the field for- 
ward to its next important step: surveying. This will open 
a new observational window onto large-scale structure and 
its evolution out to large redshifts. 



The ultimate goal of these SZ surveys is the con- 
struction of large cluster catalogs with multi-wavelength 
follow-up observations in order to perform cosmologi- 
cal studies; for example, constraining cosmological pa- 
rameters with the counts and redshift distribution (e.g., 
Barbosa et al. 1996; Eke et al. 1996; Colafrancesco et al. 
1997; Haiman et al. 2001; Holder et al. 2001; Kneissl et al. 
2001; Weller et al. 2002; Benson et al. 2002). Driving this 
effort are the particular advantages of SZ-based cluster 
catalogs (Bartlett 2000; Bartlett 2001): Firstly, SZ sur- 
veys are intrinsically efficient at finding clusters at large 
redshift, due to the surface brightness constancy of the SZ 
effect^. The thermal SZ spectrum is furthermore universal, 
the same for all clusters at any redshift^. Other emission 
mechanisms, in contrast, suffer from cosmological dim- 
ming and the need for accurate fc-corrections. Secondly, 
SZ surveys select clusters based on their thermal energy. 
Since the spectrum is the same for all clusters, the to- 
tal observable SZ flux from a cluster can be expressed in a 
frequency independent manner as the integrated Compton 
y-parameter, Y = J dVl y{n), where the integral is over 
the cluster profile (see Equation 0] below) . The y param- 
eter being the pressure integrated along the line-of-sight 
{y (X J dl nT), this then implies Y oc MT, i.e., the thermal 
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energy of the intracluster medium (ICM). This is impor- 
tant, because the total thermal energy of the ICM is given 
by energy rc-partition during cluster collapse and is inde- 
pendent of any thermal or spatial structure in the gas. It 
is hence a more robust quantity than, for example, the 
X- ray emission measure that depends in a more compli- 
cated fashion on both the ICM density and temperature. 
Hydrodynamical simulations confirm this expectation by 
showing a tight SZ flux-mass relationship with little scat- 
ter (da Silva et al. 2003). Object selection in a flux-limited 
SZ catalog is therefore relatively easy to interpret in terms 
of cluster mass and redshift. For instance, it is easy to show 
that the minimum detectable cluster mass is almost inde- 
pendent of redshift. This is particularly advantageous for 
evolutionary studies, because one is able to follow the evo- 
lution of the same kind of object over redshift, instead of 
comparing massive objects at high redshift to less massive 
ones at low redshift, as is the case with X-ray samples. 

Detailed follow-up of SZ surveys will, however, be 
time-consuming, and an enormous effort for the more than 
10^ clusters expected from Planck. Large-area photomet- 
ric surveys in the optical and infrared (e.g., SDSS) will 
help (Bartelmann & White 2002), but it is important to 
identify the kind of science that may be done directly with 
the two dimensional catalog of cluster positions and SZ 
fluxes, what we will refer to in the following as the SZ pho- 
tometric catalog. This will certainly be the first science to 
be performed. Source counts represent the primary avenue 
of 2D study that has been discussed extensively in the lit- 
erature. In this paper, we examine the next higher order 
catalog statistic, namely, the angular correlation fimction 
wie) of SZ-detected clusters (Diaferio et al. 2003). We 
quantify its information content and study its potential 
use as a cosmological probe. The angular function samples 
the catalog redshift distribution, because it is a projection 
of the real-space correlation function along the line-of- 
sight. With an appropriate model for the real-space cor- 
relation function of the catalog, we may gain some insight 
on this distribution, and hence on the underlying cosmo- 
logical model. 

Moscardini ct al. (2002) studied the 3D clustering 
properties of SZ detected clusters, taking the Planck sur- 
vey as an example and accounting for evolutionary effects 
along the past light-cone. Diaferio et al. (2003) examined 
the angular function of SZ clusters as a means of identi- 
fying probable physical cluster pairs (in 3D) and super- 
clusters. Our modeling is similar to theirs, although we 
focus instead on the angular function as a cosmological 
probe permitting the extraction of cosmological informa- 
tion from a photometric SZ catalog. The idea is of course 
not new, and has been applied in the past to, for exam- 
ple, optical and X-ray cluster catalogs; but we reiterate 
the advantages of an SZ catalog in this context: the cluster 
selection function is relatively easy to model (compared to 
other observing bands), and it extends out to large red- 
shift, giving a longer base-line for viewing evolutionary 
effects. 



We should distinguish at the outset the difference be- 
tween the angular power spectrum, Cf^, of SZ-induced 
temperature fluctuations (secondary anisotropics in the 
cosmic microwave background [CMB]) and the angular 
correlation function of detected clusters in a SZ survey, 
w{9). The angular power spectrum Cf^ is a two-point 
statistic quantifying the integrated contribution of the en- 
tire cluster population to the CMB sky brightness fluctu- 
ations. It is dominated by the Poisson term and its overall 
shape is determined by the mean SZ proflle of clusters. 
Cluster-cluster correlations add additional power on the 
order of 20 — 30% of the pure Poisson term (Komatsu & 
Kitayama 1999). Since it is deflned relative to the mean 
cosmic microwave background temperature, we expect the 
SZ fluctuation power to increase with the surface density 
of clusters on the sky. This is quantitatively confirmed 
by both numerical simulations and analytical calculations 
that indicate Cf^ oc (t| (all other factors held constant), 
where erg is the amplitude the density perturbations, the 
quantity most directly influencing cluster abundance. The 
fluctuation power spectrum is an analysis method appro- 
priate in a low signal-to-noise context (the current sit- 
uation) where individual source identification is not pos- 
sible'*. Fluctuations induced by the SZ effect have been 
invoked as a possible explanation for the excess power at 
high multipole I reported by the CBI collaboration (Bond 
et al. 2002) and consistent with new VSA data (Grainge 
et al. 2002). If this were entirely due to the cluster popula- 
tion, it would imply a surprisingly large value for crs (> 1; 
Bond et al. 2002; Komatsu & Seljak 2002; Holder 2002), 
although an important contribution from heated gas at 
reionization may also be expected (Oh et al. 2003). 

The angular function w{6) quantifies the projected 
clustering of a 2D catalog of individually detected clusters. 
It refers to the object positions and makes no reference to 
the mean background sky brightness. There arc several 
ways to imagine using the information contained in the 
SZ cluster angular function (e.g., Diaferio et al. 2003). In 
the following we choose to illustrate its use by examining 
constraints on the matter density f^M and erg in the con- 
text of flat CDM-like models. The SZ counts provide one 
constraint on a combination of these parameters. To ex- 
tract additional information from the 2D catalog using the 
angular function, we are forced to model the real-space 
cluster correlation function. In CDM scenarios, clusters 
form from peaks in the density field whose clustering may 
be analytically calculated. We adopt the approach pro- 
posed by Mo & White (1996; 2002). Any conclusions that 
we draw are, therefore, unavoidably dependent on this 
clustering model (Moscardini et al. 2002; Diaferio et al. 
2003); however, it is well founded in the context of CDM 
cosmogonies and compares well with the results of nu- 
merical simulations, at least at redshifts lower than ~ 10 



* Either for low signal-to-noise observations or when push- 
ing constraints on source counts below the detection threshold. 
It is the SZ equivalent of the "background fluctuation analysis" 
of radio and X-ray Astronomy. 
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(Reed et al. 2003; Jenkins et al. 2001). Another important 
issue concerns the modeling of the intracluster medium. 
Moscardini et al. (2002), for example, find non-negligible 
dependence of the 3D correlation function on ICM prop- 
erties. We discuss the issue below for the angular function, 
but we note here that it tends to be dominated by massive 
clusters (depending of course on the exact flux limit; see 
below) that follow relatively well observed scaling laws. 
As we show, adding constraints from the local X-ray clus- 
ter abundance permits us to simultaneously constrain the 
cosmological parameters and the cluster baryon content. 

In the following section, we give our master equation 
for calculating 'w{9) and identify the necessary modeling 
ingredients. We outline our cluster model in Section 3. 
Results for the angular function are presented in Section 
4, where we apply the results to constraining the cosmo- 
logical parameters cts and Qm- We discuss the influence 
of ICM physics on the results and how to use additional 
information from the local cluster abundance to simul- 
taneously constrain this physics and the cosmology. Our 
fiducial example is the Planck mission. Section 5 closes 
with a final discussion and summary. 

2. The Angular Correlation Function 

In this section we relate the angular correlation function 
to the real-space correlation function in the context of SZ 
observations. The 3-dimensional (auto)correlation func- 
tion ^ quantifies the 2-point clustering of a population in 
terms of the probability in excess of Poisson of finding two 
objects at a separation r. The angular correlation w{9) of 
the same population is the projection of ^ onto the sky: 



dnidn2 / dzidz] 



dnidn2^^[l 
dVi 



- HO)] 
dV2 



(1) 



dMidM2 



dflidzi dVl2dz2 
dni dn2 



dMi dM2 



(1 + C) 



where the integrals concern two lines-of-sight (los) of 
solid angles dili and dfl2 separated on the sky by an- 
gle 9. In this expression, dn/dM represents the cluster 
mass function, and we assume the small angle approxima- 
tion here and throughout. The surface density of sources 
(the counts) with integrated Compton parameter larger 
than Yiini (the SZ 'flux' limit; see the following section. 
Equation 2)) is given by 



S(i1im) 



dz 



dV 
dQdz 



M„„(lii„,z) 



flT) 

dM — (M,z) (2) 
dM^ ' ' ^ ' 



The correlation function ^ depends on cluster mass, red- 
shift and, according to statistical isotropy, on physical 
separation r{9) — r'l{zi) + ^'2(2:2) — 2rir2 cos(6'), where 
ri and r2 are angular-diameter distances. Assuming that 
correlations fall off sufficiently rapidly with distance, as 
is observed, we may take the two clusters to be at ap- 
proximately the same redshift and write ^[Mi, M2, z, r{9)]. 
We furthermore adopt a linear biasing scheme in 



which £_{Mi,M2,z,r) = b{Mi, z)b{M2, z)£_dmiz,r), where 
^dm{z,r) is the correlation function of the underlying 
cold dark matter and b{M, z) is the bias factor for 
clusters (see below). Then, using the short-hand no- 
tation 0(M,z, YiinO ^ {l/T.)b{M,z){dV/dzdn)dn/dM for 
the joint distribution of clusters in mass and redshift, 
weighted by the bias factor, we arrive at the expression 

w{9,Y,^) = 



dz I I dMidM2 (j){Mi, z)(j){M2, z) y. 

J JMum(Ynm,z) 

dz 

dr—{z)^dm[z,ri9)] 
dr 



j dz$^(z,Yiim) j dr—{z)(,Am[z.,': 



a sort of Limber's equation appropriate for SZ sources 
(Limber 1953; Peebles 1993; Diaferio et al. 2003) in which 
we explicitly show the dependence on limiting flux Yw-^. 
From this equation we clearly see that the three key ingre- 
dients are the the mass-limit function Miim(z, yjim), the 
distribution function (j) and the correlation function ^dm- 
We now discuss our modeling of each. 

3. The SZ Population 

The SZ cluster population inherits its properties from two 
sources: its constituent dark matter halos, whose proper- 
ties are the sole result of gravitational evolution, and the 
relationship between observable SZ flux and these halos, 
governed by more difficult to model baryonic physics. It 
is reasonable to characterize dark halos by their mass and 
redshift, and we will apply the results of N-body simula- 
tions that give both their abundance (i.e., the mass func- 
tion) and spatial correlations (i.e., clustering bias b{M,z) 
and the correlation function ^dm) as a function of these 
two fundamental descriptors. More difficult to model, the 
baryonic physics of the cluster gas requires particular at- 
tention to various (and at times contradictory) observa- 
tional constraints and theoretical scaling laws. 

3.1. Halo properties 

The abundance of galaxy clusters is given by the mass 
function of collapsed objects, which is completely speci- 
fied once the linear power spectrum of dark matter per- 
turbations is specified. For the latter we adopt the BBKS 
(Bardeen et al. 1986) transfer function (see also below), 
while for the mass function we employ the fitting formula 
(improved Press-Schechter) given by Sheth & Tormen 
(1999): 



— (M, z)dM = A^^^)hiim'''A — 



dMZ) 



where p is the universal mean mass density and the con- 
stants A « 0.322 and q = 0.3; the parameter v' = -Jav, 
where v = ntJl/iK.r\ is the usual critical peak height 



D{z)a{M) 
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(Sc ~ 1.69) normalized to the mass-density perturbation 
variance a{M) in spheres containing mass M, and the con- 
stant a = 0.707. The expression for the growth factor for 
flat models with A > is taken from Carroll et al. (1992): 



D{z, Q,m,^Ia) = 



a(o){i+z) 

-(l + r!M(2)/2)(l + r!A(z)/70)]- 



with the definitions fluiz) = flui^ + /E'^{z), ^a{z) = 
QA/E^iz), and E^{z) = [f^A + (1 - ^^m - ^A)il + zf + 
^^m(1 + zf'\ \ r^A/ and J^a written without an explicit red- 
shift dependence will indicate present-day values (z = 0). 

As mentioned above, we use a linear bias scheme to re- 
late the cluster-cluster correlation function to that of the 
dark matter and employ the analytic fitting formula for 
6(M, z) given by Sheth et al. (2001); the formula includes 
corrections for ellipsoidal perturbation collapse: 



6(M, z) = 1 + 



1 



^fa&ciz) 



[Va(c 



\fah(av 



2\l-cl 



{av^Y + h{l-c){l~c/2) 

where 5c, v and a are given above, and h = 0.5 and c — 0.6. 
We model the linear dark matter perturbation spectrum 
with the BBKS transfer function with shape parameter 
fixed atV = 0.25 and scale-invariant primordial density 
fluctuations in — 1). This seems to provide a good fit to 
galaxy clustering data (Percival et al. 2001) and is consis- 
tent with constraints on n from CMB anisotropies (Spergel 
et al. 2003, and references therein). The resulting linear 
theory ^dm is adequate on most scales (Q > 10'), although 
we also include non-linear corrections according to the 
fitting formula developed by Peacock & Dodds (1996). 

3.2. Intracluster Medium 

With the abundance and clustering of halos now specified, 
we next relate the observable SZ flux to cluster mass and 
redshift. This relation is particularly robust from a theo- 
retical viewpoint, contrary to, for example, the situation 
for X-ray luminosity. In the non-relativistic regime, the 
surface brightness of the thermal SZ effect - measured rel- 
ative to the mean sky intensity ~ at position f2 on a cluster 
image is 

where the Compton y parameter is an integral of the pres- 
sure along the line-of-sight 



dl 



kT 



los h mc^ 



jnar 



with T the gas temperature (strictly speaking, that of 
the electrons), ctt the Thompson cross section, k and m 
the Boltzmann constant and the electron mass, respec- 
tively, and where j^, is a universal spectral function that 



SZ Angular Correlotion Function 
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Fig. 1. Angular correlation function of SZ detected clus- 
ters calculated in two different flat models (JIm = 0.3 and 
(78 — 0.9 as the dashed hnes, and ilu = 1 and as — 0.7 
as the solid lines) for three flux limits (liim — 10~^, lO^'', 
10~^ arcmin^, decreasing from top to bottom). The mod- 
els adopt the same linear matter power spectrum (shape 



parameter F = 
(see text). 



0.25) and = 7.6 x IQ-^h"^^^ arcmin^ 



is the same for all clusters, independent of their proper- 
ties (Birkinshaw 1999). Since the thermal SZ spectrum is 
the same for all clusters, we may express the total flux 
in a frequency independent manner using the integrated 
Compton parameter Y — J dfl y(ri), where the integral is 
over the cluster profile. The total flux density (e.g., in Jy) 
is then Si, = Y j^. This is the total flux (density) measured 
by an experiment with low angular resolution in which 
clusters are simply point sources. From these definitions, 
we find that the observable flux is directly proportional to 
the total thermal energy of the ICM (e.g., Barbosa et al. 
1996): 



y(M, z) 



k(jT N^T 



fg^{M,z)T{M,z)M 



(4) 



where is the total number of electrons, /gas is the ICM 
mass fraction, Z?ang('Z) = Ho^'^angiz) is the angular di- 
ameter distance, and T is to be understood as the mean 
(particle, and not emission-weighted) electron tempera- 
ture; note that with this understanding, the relation does 
not depend on any assumption of isothcrmality. It is this 
direct relation between observable SZ flux and thermal gas 
energy that lies at the heart of some of the advantages of 
SZ over X-ray surveys (Bartlett 2001). 
We have taken care to write the quantities /gas and T 
as general functions of mass and redshift. In the absence 
of efficient heating/cooling and gas reprocessing, the clus- 
ter population will be fully self-similar. Simple theoretical 
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arguments based on energetics during collapse suggest in 
this case the existence of a scaling law between cluster 
temperature and mass: 



1 



A(z) 



(5) 



adopting the notation of Pierpaoli et al. (2002). In this ex- 
pression, the mass Afis is measured in units of 10^^ M©, 
A(2;) is the full non-linear overdensity inside the virial 
radius relative to the critical density 178). In the self- 
similar model, the gas mass fraction /gas is constant, es- 
sentially proportional to the universal ratio rJe/^^M- 
These scaling expectations are indeed manifest in hydro- 
dynamical simulations that neglect cooling, and supported 
by observations of the more massive clusters with cooling 
timescales longer than the Hubble time (Mohr et al. 1999; 
Allen et al. 2001; Finogucnov et al. 2001; Xu et al. 2001; 
Arnaud et al. 2002). Although there is good indication 
that clusters are a population with rather regular proper- 
ties, as suggested by these scaling laws, there is also direct 
evidence that it is not exactly self-similar. Deviations are 
most pronounced for the lower mass objects (T < 2 keV), 
as perhaps expected since they have shorter cooling times 
and plausible energy injection mechanisms more readily 
compete with their gravitational energy (Ponman et al. 
1999; Lloyd-Davies et al. 2000). More generally, one may 
write 



Y{M,z) = Y,,{z)Mli^+'^{l + zr 



(6) 



The factor Yi^{z) incorporates all the z dependence of 
Eqs. H4I5I) and is defined such that the self-similar model 
corresponds to a = 7 = 0. 

Although we have discussed Eq. © via scaling rela- 
tions for T and /gas, it is more pertinent to consider it as 
a single relation between observable SZ flux and cluster 
mass and redshift, one that expresses the proportional- 
ity between the observable flux and the gas thermal en- 
ergy. Hydrodynamical simulations incorporating cooling 
and pre-heating (da Silva et al. 2003) indicate a between 
0.1 and 0.2 and 7 ~ 0, with a very tight scatter about the 
relation, much tighter than corresponding relations for X- 
ray quantities. This confirms our expectations that the SZ 
flux is a more robust quantity than its X-ray counterparts. 
Even though very low mass clusters are included in the fit, 
only mild deviations from self-similarity in this relation 
are seen in the simulations. These deviations are most pro- 
nounced in the low mass systems, while the higher mass 
systems appear compatible with the self-similar scaling 
laws for T-M and /gas- Caution is still warranted, how- 
ever, due to the unstable nature of cooling, issues of nu- 
merical resolution, and the fact that these codes model 
the gas as a single phase medium; the simulation volume 
also contains mostly low mass systems, with only a few 
clusters with A/15 > 0.2. Simulation results are nonethe- 
less roughly consistent with X-ray observations (Borgani 
et al. 2002; Muanwong et al. 2002). 

At our fiducial Planck flux limit oiY = 10"** arcmin^, 
we have checked that the counts and value of the angular 



function at 30 arcmins are little affected by clusters below 
T = 2 keV, at least at ^Im < 0.6. We therefore concen- 
trate on the self-similar case with a = 7 = 0. On the 
other hand, we consider the normalization of Eq. as 
a parameter free to vary within certain limits suggested 
by X-ray observations of /gas and T. Moscardini et al. 
(2002) have argued that this freedom makes it difficult 
to use only SZ observations to constrain cts and Qm- To 
overcome this modeling uncertainty, we combine SZ obser- 
vations of both the number counts and the angular cor- 
relation function with constraints arising from the local 
abundance of X-ray clusters. We shall find that the three 
kinds of observations are complementary and lead to con- 
straints on the cosmological parameters. 
To give a feel for the order of magnitude, we note that 



Yi5{z) = (7.4 X 10-^/i^/s 



( 



178 



arcmm 



1/3 r 



( kcV ) \, 0.07h--V2 



178 



1 

1/3 



) Oa(z) 



1 



o S^a(z) 

^ A(2) 



(7) 



For reference, T, — 1.2 according to the simulations of 
Evrard et al. (1996) and /gas = Omh'^-^ from Mohr et 
al. (1999). At an observation frequency of 2.1 mm, the 
maximum decrement of the thermal spectrum, slY ~ 7.4 x 
10~^ arcmin^ corresponds to a flux density of ^ 7 mJy. In 
all the analysis a minimum cluster mass of lO^^h^^ Mq is 
imposed. 

4. Results 

4.1. The angular correlation function 

Calculated angular correlation functions are shown in 
Figure n for three different limiting flux values (ilim) in 
two different flat cosmologies {Hm = 0.3 and JIa — 0.7, 
with (Tg ~ 0.9, and flu = 1 with erg = 0.7; e.g., 
Blanchard et al. 2000); the Y~M normalization is = 
7.6 X IQ-^h'^/^ arcmin . Our angular correlation function 
is consistent with the results of Diaferio et al. (2003) at 
a separation of 30 arcmins. As shown by the latter au- 
thors, small scales are affected by the specific choice of 
clustering bias function b{M, z), but the model differences 
drop to ^ 10 % at 30 arcmins. This is comparable to the 
statistical measurement errors expected in the case of the 
Planck survey, as discussed below. We therefore use the 
angular function at 30 arcmin separation in our analysis 
of cosmological parameters. 

Figure |2] shows the cluster selection function $(2;), de- 
fined in Eq. (|3J), at — 10^* arcmin^ in the two cos- 
mologies and for various normalizations erg. The broad 
break in the angular function in Figure ^ corresponds to 
the break in ^dm just beyond ~ 10/i~^ Mpc. For exam- 
ple, the break occurs ~ Ideg for Yjim = 10^"* arcmin^ 
in the critical model. According to Figure [21 the selection 
function peaks around z ~ 0.2 (for the chosen value of 
as — 0.7), projecting the break to an angular separation 
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Fig. 2. The selection function <i> of SZ-detected clusters with Y > Yjim = 10^^ arcniin^ (see Eq. O as a function of 
redshift z and for different values of erg in the low-density model (J^m ~ 0.3, on the left) and in the critical model 
(r^M = 1, on the right). From top to bottom, the curves are for increasing as — 0.5, 0.6, 0.7, 0.8, 0.9, 1.1, 1.4. The models 
adopt the same linear matter power spectrum (shape parameter T = 0.25) and Y* — 7.6 x lO^^h'^^^ arcmin^. 



of 6* ~ 10/i-V(3000/i-i X 0.2) - Ideg. Figures d and H 
visually illustrate how the angular function encodes infor- 
mation on the catalog radial distribution. 

Although the two models have comparable angular 
correlations at the bright end (the upper curves in Fig. ^ 
corresponding to Y^^ = 10^^ arcmin^), their dependence 
on catalog depth clearly differs. At the bright end, 
we are mainly observing the local cluster population, 
which is essentially the same in both models since the 
present-day abundance is the same and the density 
perturbation power spectrum is fixed (F = 0.25). Note 
that the low-density model has a slightly steeper slope at 
small separation, due to its greater non-linear evolution. 
The angular function of the low-density model decreases 
and shifts to the left more rapidly with survey depth 
than in the critical model. The overall trend is easily 
understandable and due to the fact that the selection 
function broadens and peaks at higher redshift with 
survey depth, moving correlations to smaller angular 
scales and generally washing out the signal as more 
clusters are projected along the line-of-sight. In bright 
galaxy surveys, which sample the local universe where 
space is approximately Euclidean and galaxy evolution 



may be ignored, the dependence of the angular function 
on survey depth follows an important scaling law that is 
independent of the underlying cosmological model. No 
such universal scaling law obtains in the SZ case, because 
the radial distribution extends out to large redshifts. For 
the relevant flux limits, a SZ catalog therefore samples 
evolution in the cluster population. Since the cluster 
population evolves less rapidly with redshift in the low- 
density model, the angular correlation function therefore 
shifts down more rapidly with survey depth than in 
the case of the critical model. The angular correlation 
function of SZ clusters is therefore a cosmological probe. 



4.2. Combining the SZ angular function and counts 

As an illustration of this probe, we next consider the 
dependence of the angular function and of the SZ counts 
on the cosmological parameters {ilM,<^8)- In Figure |31 
we show the predicted counts and angular correlation 
function for a set of flat models over a grid in the 
(cTg, riM)~plane and for different limiting flux values, 
liim = 10~^ arcmin^ (top left), Yu^ 
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(top right), and Fiim = 10~^ arcmin^ (bottom left). 
Contours of the angular correlation function are shown as 
continuous lines and refer to its value at a separation of 
9 = 30 arcmins. The counts contours at the correspond- 
ing flux limits are given as dotted lines and labeled in 
units of 1/sq.deg. This figure is constructed for a self- 
similar cluster population with a Y-M normalization, 
K = 7.6 X lO-^/i'^/^ arcmin^ (see Eq. [7|); the effects of 
varying ICM properties will be discussed below. 
We see from the figure that the counts and angular 
function carry more complementary information as the 
catalog increases in depth. At the high flux limit of 
^lim = 10~^ arcmin^, the catalog primarily probes the 
local universe; the two sets of contours are parallel and 
there is no information permitting us to constrain the two 
cosmological parameters. As we move deeper, the two sets 
of contours begin to cross, indicating the we are obtaining 
useful cosmological information capable of constraining 
the two parameters. The limit Yuai = 10~^ arcmin^ 
is representative of the Planck mission (Aghanim et 
al. 1997; Bartelmann 2001; and references therein), 
while Yiim = 10~^ arcmin^ would correspond to deeper 
ground-based surveys (probably attaining the confusion 
limit). 



gular correlation function (Eq. |S| below). We therefore fo- 
cus on the effects of changing the normalization, Eq. Q, 
of the self-similar SZ flux-Mass relation (a = 7 = 0). 
Physically, this represents an uncertainty in the average 
thermal gas energy of galaxy clusters. Figure 0] shows the 
effect of changing at a flux limit of Yjim = lO""' arcmin^. 
From the figure, we see that changing this normalization 
stretches the counts and angular function contours verti- 
cally by roughly the same amount, their point of intersec- 
tion (if any) remaining at roughly the same value of ^m- 
For example, a factor 2 change in normalization moves 
the intersection by ^ 10 — 20 % in as (e.g., at Qm ^ 0.3). 
This represents an inherent systematic caused by mod- 
eling uncertainty associated with the ICM. Uncertainties 
related to ICM modeling have been extensively discussed 
by Moscardini et al. (2002) for the full 3D spatial correla- 
tion function of SZ clusters. 

Including the uncertain SZ flux normalization, Y^, we now 
have three parameters to determine, and so additional in- 
formation is needed; for example, an observational pro- 
gram to determine the Y-M relation on a representative 
sample of clusters. Another tactic is to use the constraint 
arising from the local cluster abundance as measured by 
the X-ray temperature function. 



4.3. Influence of the ICM 

The single relation Eq. 10) incorporates all of the ICM 
physics, and from the nature of the SZ fiux, we expect 
it to be rather tight with little scatter. There is virtually 
no direct observational information on this relation, while 
numerical simulations (da Silva et al. 2003) do indeed con- 
firm a small scatter. They also indicate that there is only a 
slight deviation from the self-similarity with a = 0.1 — 0.2 
and 7 = down to masses well below 10^'' M©. X-ray 
observations, on the other hand, demonstrate deviations 
from self-similarity. This is most clear for low-mass sys- 
tems with temperatures below T ^ 2 keV, but even in 
the richer systems the observed L-T deviates from the 
self-similar expectation. These results are beginning to 
be reproduced by numerical simulations including cool- 
ing/heating (e.g., Borgani et al. 2002), including those 
that indicate only slight deviation from self-similarity for 
the SZ relation Eq. ® (Muanwong et al. 2002). That X- 
ray observables show greater deviation from self-similarity 
than the SZ flux is not surprising, given that the former 
are more sensitive to spatial/temperature structure in the 
gas than the latter, as emphasized previously. This is pre- 
cisely the advantage of SZ surveys. 

As mentioned, low mass systems with the greatest 
deviation from self-similarity in X-rays contribute lit- 
tle to counts and the angular function at a flux limit of 
^iim = 10"'' arcmin^. For example, for Hm < 0.6 and val- 
ues of (Tg in the range 0.6-1, the uncertainty on as due 
to the presence of low mass systems with T < 2 keV is 
around 5%, well within our estimated errors on the an- 



4.3.1. Adding the local cluster abundance constraint 

For a flxed value of , the local X-ray temperature func- 
tion constrains the cosmological parameters erg and S7m 
to a well-defined curve in the plane: cts^^m^ ~ 0.6r^"''-^ 
(Pierpaoli et al. 2002). The exact relation depends some- 
what on the chosen mass function. Consider such a con- 
straint, shown as the middle dashed line in Figure 
Supposing that we measure now both the SZ counts 
and angular function, finding 7.7 sq. deg and w{0 = 
30 arcmins) = 0.024, respectively, ai Y = 10""* arcmin^ 
and given as the solid contours in the figure. Apart from 
the cosmological parameters as and ^m, the SZ contours 
depend only on K,, whose value is uncertain; by changing 
y*, we move the SZ contours around the plane, as previ- 
ously described. However, there is a unique value for Y^ 
that reduces the intersection of the three lines to a sin- 
gle point in the plane. By adjusting Y^, to this value, we 
simultaneously constrain the cosmological parameters and 
the normalization Y^ . In the figure, this intersection lies on 
top of the true underlying model. The determined value 
of together with the known imply a constraint on 
/gas^- By adding the information on the local cluster abun- 
dance, we have pinned down the relevant ICM physics and 
constrained the cosmological parameters, eliminating the 
primary uncertainty in using SZ clustering observations 
(Moscardini et al. 2002). Unfortunately, T, is at present 
only poorly known, with simulations and observations in- 

^ We are implicitly assuming that mean electron tempera- 
ture, relevant to the SZ effect, and the observable X-ray tem- 
perature, parametrized by T*, are the same; this need not be 
the case. 
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Fig. 3. Contours of the SZ angular correlation function arc shown as solid lines in the {as, riM)~planc for different lim- 
iting flux values: top left, Yjim = 10~^ arcmin^; top right, Yjim — 10~* arcmin^; and bottom left, Yu^ — 10~^ arcmin^. 
The Hubble constant is fixed at /i = 0.7. Each contour is labeled with the value of the correlation function at a 
separation of = 30 arcmins. Contours of the SZ cluster counts at the same flux limits are shown as the dotted lines, 
and labeled in units of 1/sq.deg. The SZ flux normalization = 7.6 x lO^^h'^/^ arcmin^ . 



dicating values generally in the range 1.2 — 2 (Pierpaoli et 
al. 2002; Muanwong et al. 2002; Huterer & White 2002). 
We expect that it will be much better determined by the 
time large SZ catalogs become available, making this kind 
of analysis possible. 

This example illustrates the utility of the angular 
function: with only the local cluster abundance and 
the SZ counts, we cannot determine the cosmological 



parameters erg and ^Im, due to the uncertainty on y*. 
Adding the angular function breaks the degeneracy. We 
are thus able to constrain the cosmology with the 2D 
SZ catalog, without recourse to redshift determinations. 
Constraints obtained in this way are particularly useful 
for their complementarity to constraints from CMB 
anisotropy and distance measurements with SNIa. For 
example, with SZ clusters we measure as directly on 



S. Mei and J.G. Bartlett: Cosmological information from a 2D SZ catalog 



9 




1.10 F 




0.50 F , , , ''I. , , , I , I , , , I , , V-^ 
0.0 0.2 0.4 0.6 0.8 1.0 

n 



Fig. 4. Effect of varying F*. Tlie contours liave the same meaning as in tlie previous figure, and h = 0.7 as before. All 
three panels correspond to a flux limit of Yni^ — 10"* arcmin^, but changes from 3.8 x lO^^/i^/® arcmin^ (upper left) 
to 7.6 X 10~^/i.^/^ arcmin^ (upper right; same as upper right panel of previous figure) and finally 1.5 x 10~*/i^/® arcmin^ 
(lower left). The contours and their spacing are stretched upward as Y^, decreases. The intersection points of the two 
sets of contours essentially move vertically, changing in erg and remaining roughly fixed in J^m ■ 



the relevant scales and constrain Qm, as opposed to the 
physical density iluh'^ in the case of the CMB. 

As mentioned above, clusters below ~ 2 keV con- 
tribute little to the contours at our fiducial flux limit of 
Y = 10^"* arcmin^. At lower flux limits, more representa- 
tive of future ground-based surveys, we could expect that 
the possible effects of deviations from self-similar scaling 



become more important, although simulations at present 
indicate only mild effects on the SZ flux. 

4.4. Discussion of statistical errors 

Before concluding, we briefly examine the statistical er- 
rors expected on a measurement of the SZ angular cor- 
relation function. Although a thorough analysis of this 
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Fig. 5. Example of constraints in flat models. The true 
underlying model corresponds to flu = 0.3, cs = 0.9, 
/gas = and T, — 1.5 keV, corresponding to 

y* = 1.1 X lO^^h"^^^ arcmin^. For this model, one would 
observe 7.7 clusters/sq.deg. for the counts and w — 0.024 
for the angular function at 30 arcmins, respectively, both 
at a flux limit oiY ~ 10"'' arcmin^. The bold (heavy blue) 
line shows the constant counts contour in the plane, while 
the thin (lighter blue) line indicates the angular function 
contour at the observed value. The middle dashed line rep- 
resents the constraint from the local X-ray cluster tem- 
perature function, for the true value of T*; the other two 
dashed lines correspond to the constraints for = 2 keV 
(bottom dashed curve) and T* = 1 keV (upper dashed 
curve). The three sets of contours cross at a unique point. 
The gray shaded band indicates the estimated statistical 
uncertainty on w attainable with the Planck survey (the 
error on the counts is roughly the thickness of the con- 
tour) . 



issue requires detailed simulations of a given survey, we 
may nonetheless make some general arguments to gain 
insight into what may eventually be achieved. We take 
the Planck survey as our example. The resolution of the 
Planck SZ catalog will be on the order of 5 arcmins (at 
best), so we can only expect to measure the correlations 
on larger scales (as the 30' we have chosen in our analysis), 
and the flducial sensitivity expected is F ~ 10""* arcmin^ 
(Aghanim et al. 1997; Bartelmann 2001). Since the an- 
gular correlations are small in this context (< 0.1), we 
may estimate the (statistical) error on a measurement of 
w^O) as the Poisson variance in the number of pairs, ripair, 
at this separation (e.g., Peebles 1980; Landy & Szalay 
1993). This quantity is determined by the counts as fol- 
lows: Suppose that we measure w in a annular bin of 
width A9 at angular distance 9 from a cluster. The mean 



number of clusters in this ring is (n) — 27r0A6'E(yiim), 
from which we deduce that the total number of pairs 
at this separation in a catalog of N clusters is about 
"pair « (1/2) X (n) = Ntt9A9T,. In other words, we 
estimate the error to be 



Awi9,Ynn.) 



~ ^ P(^lim)] ^ I (4^ 

2.78 X 10-3 (iJgs) (M 



-1/2 



1/2 / s 

ydcg^ 



(8) 



Notice that this statistical error depends essentially on the 
number counts S. This leads to the gray shaded area in 
Figure |31 

5. Discussion and Conclusion 

We have calculated the angular correlation function of SZ- 
detected clusters in order to evaluate its usefulness for ex- 
tracting cosmological information directly from a 2D SZ 
cluster survey, before 3D follow-up. The angular correla- 
tion function of SZ detected clusters differs from the an- 
gular power spectrum of SZ induced CMB anisotropics, 
which is dominated by the Poisson term. The different 
scaling of angular correlations with survey depth visu- 
ally demonstrates the cosmological sensitivity of the an- 
gular function (See FigsnJID). illustration, we consid- 
ered the parameter pair (r^MiCg) in the context of flat 
A-CDM models. We found that at sufficient depth (e.g., 
llim ^ lO""* arcmin^, comparable to the Planck mission), 
the counts and angular function combined can constrain 
these parameters. Modeling uncertainty associated with 
the ICM may be reduced by measuring the T-M rela- 
tion and adding the corresponding constraint from the 
local abundance of X-ray clusters; in this way, the two 
cosmological parameters and the SZ flux normalization 
y* may be found. Deeper ground-based surveys (e.g., 
yiim ~ 10"^ arcmin^), will pick up a larger number of low- 
mass objects whose ICM properties require more careful 
modeling of deviations from self-similarity. 

The accuracy and precision with which one will be 
able to measure the angular function (and the counts) is 
clearly an important issue. We only briefly touched on this 
point with a simple estimate of the expected statistical er- 
rors expected on ■w{9) in the case of the Planck mission. 
A more detailed examination incorporating simulations of 
the SZ survey characteristics, including instrument noise 
and foreground contamination, as well as of the catalog ex- 
traction algorithms, is needed. The survey selection func- 
tion will in reality depend on these details (Bartlett 2001; 
Mehn et al. 2003; White 2003), and so wiU the errors on 
the counts and measured angular function. 

It is clear that the most powerful constraints from an 
SZ survey will come from its measured redshift distribu- 
tion. The interesting point is, however, that the method 
proposed here increases the immediate scientific return 
from an SZ survey by offering a way to obtain pertinent 
cosmological constraints using only a 2D SZ catalog, with- 
out recourse to follow-up observations. 



S. Mei and J.G. Bartlett: Cosmological information from a 2D SZ catalog 



11 



Acknowledgements. We are grateful to Lauro Moscardini, 
Antonaldo Diaferio and Antonio da Silva, and the anonymous 
referee for their very useful comments that helped the improve- 
ment of this paper. We would also like to thank Antonella De 
Luca and Milan Maksimovic, who first opened the contacts 
between us. S. Mei acknowledges support from the European 
Space Agency External Fellowship programme. 



References 

Aghanim, N., de Luca, A., Bouchet, F.R. et at 1997, A&A, 
325, 9 

Aghanim, N., Castro, P.G., Melchiorri, A. et al. 2002a, A&A, 
393, 381 

Aghanim, N., Hansen S., Pastor S. et al. 2002b, 

astro-ph/0212392 
Allen, S. W., Schmidt, R.W., Fabian A.C., 2001, MNRAS, 328, 

L37 

Arnaud, M., Aghanim, N., Neumann, D. M. 2002, 389, 1 
Barbosa D., Bartlett J.G., Blanchard A. et al, 1996, A&A, 
314, 13 

Bardeen, J.M., Bond, J.R., Kaiser, N., Szalay, A.S. (BBKS) 

1986, ApJ, 304,15 
Bartelmann, M. 2001, A&A, 370, 754 
Bartelmann, M. & White, S.D.M. 2002, A&A, 388, 732 
Bartlett J.G. 2000, astro-ph/0001267 

Bartlett J.G. 2001, review in "Tracing cosmic evolution with 
galaxy clusters" (Sesto Pusteria 3-6 July 2001), ASP 
Conference Series in press |astro-ph/01112irj 

Benson, A.J., Reichardt, C, Kamionkowski, M. 2002, MNRAS, 
331, 71 

Benson, A.J., Church, S.E., Ade, P. A. R. et al. 2003, 

astro-ph/0303510 
Birkinshaw M. 1999, Proc. 3K Cosmology, 476, American 

Institute of Physics, Woodbury, p. 298 
Blanchard A., Sadat R., Bartlett J.G. et al. 2000, A&A 362, 

809 

Bond J.R., Contaldi C.R., Pen, U-L. et al. 2002, submitted to 

ApJ, astro-ph/0205386 
Borgani S., Governato F., Wadsley et al. 2002, MNRAS 336, 

409 

Carlstrom, J.E., Holder, G.P., Reese, E.D. 2002, ARA&A, 40, 
643 

Carroll, S.M.; Press, W.H., Turner, E. L. 1992, ARA&A, 30, 
499 

Colafrancesco, S., Mazzotta, P., Vittorio, N. 1997, ApJ, 488, 
566 

da Silva, et al. 2003, in preparation 

Diaferio, A., Nusser, A., Yoshida, N. et al. 2003, MNRAS, 338, 
433 

Eke, V.R., Cole, S., Frenk, C.S. 1996, MNRAS, 282, 263 
Evrard, A.E., Metzler, C.A. & Navarro, J.F. 1996, ApJ 469, 
494 

Finoguenov, A, Reiprich, T.H. & Bohringer, H. 2001, A&A 
368, 749 

Grainge K., Carreira, P., Cleary K. et al. 2002, 

astro-ph/0212495 
Grego, L., Carlstrom, J.E., Reese, E.D. et al. 2002, ApJ, 552, 

2 

Haiman, Z., Mohr, J. J., Holder, G. 2001, ApJ, 553, 545 
Holder, G., Haiman, Z., Mohr, J.J. 2001, ApJ, 560, Llll 
Holder G.P. 2002,ApJ, 578, LI 



Holzapfel, W.L., Ade, P.A.R., Church, S.E. et al. 1997, ApJ, 
481, 35 

Huterer, D. & White, M. 2002, ApJL, 578, L95 
Jenkins A., Frenk C.S., White S.D.M. et al. 2001, MNRAS, 
321, 372 

Jones M.E., Edge A.C., Grainge, K. et al. 2001, 

astro-ph/0103046 
Kneissl, R., Jones, M.E., Saunders, R. et al. 2001, MNRAS, 

328, 783 

Komatsu, E., Kitayama, T. 1999, ApJ, 526, LI 
Komatsu, E. & Seljak, U. 2002, MNRAS, 336, 1256 
Landy, S.D. & Szalay, A.S. 1993, ApJ, 412, 64 
Limber D.N. 1953, ApJ 117, 134 

Lloyd-Davis E.J., Ponman T.J. & Cannon D.B. 2000, MNRAS 
315, 689 

Mason B.S., Myers S.T. & Readhead A. C.S. 2001, ApJ 555, 
Lll 

Melin, J.-B. et al. 2003, in preparation 

Mo, H.J., White, S.D.M. 1996, MNRAS, 282, 347 

Mo, H.J., White, S.D.M. 2002, MNRAS, 336, 112 

Mohr, J.J., Mathiesen, B., Evrard, A. E. 1999, ApJ 517, 627 

Moscardini, L., Bartelmann, M., Matarrese, S. et al. 2002, 

MNRAS, 340, 102 
Muanwong, O., Thomas P.A., Kay, S.T. et al. 2002, MNRAS, 

336, 527 

Oh, S.P., Cooray, A., Kamionkowski, M. 2003, 

astro-ph/0303007 
Peacock J. A., Dodds, S. J. 1996, MNRAS, 280, L19 
Peebles P.J.E. 1980, The Large-Scale Structure of the 

Universe, Princeton Series in Physics, Princeton University 

Press (Princeton: New Jersey) 
Peebles P.J.E. 1993, Principles of Physical Cosmology, 

Princeton Series in Physics, Princeton University Press 

(Princeton: New Jersey) 
Percival W.J., Baugh CM., Bland-Hawthorn J. et al. 2001, 

MNRAS 327, 1297 
Pierpaoh E., Borgani S., Scott D. et al. 2002, astro-ph/0210567 
Ponman T.J., Cannon D.B., Navarro J.F. 1999, Nature 397, 

135 

Reed, D., Gardner J., Quinn T. et al. 2003, astro-ph/0 301270^ 
Reese, E.D., Carlstrom, J.E., Joy, M. et al. 2002,ApJ, 581, 53 
Sheth R.K. & Tormen G. 1999, MNRAS, 308, 199 
Sheth R.K., Mo H.J. & Tormen G., 2001, MNRAS, 323, 1 
Spergel, D.N., Verde L., Peiris H.V. et al. 2003, 

astro-ph/0302209 
Sunyaev, R. A. & Zel'dovich, Ya. B. 1970, Comments 

Astrophys. Space Phys., 2, 66 
Sunyaev, R. A. & Zel'dovich, Ya. B. 1972, Comments 

Astrophys. Space Phys., 4, 173 
Weller, J., Battye, R.A., Kneissl, R. 2002, PhRvL, 88, 1301 
White M. 2003, astro-ph/0302371 
Xu, H., Jin, G.,Wu, X.-P. 2001, ApJ 553, 78 



